Stabilizing topological phases in graphene via random adsorption 
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We study the possibihty of reaUzing topological phases in graphene with randomly distributed ad- 
sorbates. When graphene is subjected to periodically distributed adatoms, the enhanced spin-orbit 
couplings can result in various topological phases. However, at certain adatom coverages, the inter- 
valley scattering renders the system a trivial insulator. By employing a finite-size scaling approach 
and Landauer-Biittiker formula, we show that the randomization of adatom distribution greatly 
weakens the intervalley scattering, but plays a negligible role in spin-orbit couplings. Consequently, 
such a randomization turns graphene from a trivial insulator into a topological state. 
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Introduction — Due to its unusual structure, graphene 
displays two other internal degrees of freedom in addition 
to the real spin: AB sublattice and valley [1 . The for- 
mer refers to the interpenetrating triangular sublattices 
A and B, while the latter refers to a pair of inequiva- 
lent Dirac cones K and K' in the Brillouin zone. When 
graphene is deposited on substrates or adsorbed with 
heavy atoms, the interaction-induced symmetry break- 
ing can open bulk energy gaps to support a rich va- 
riety of topological phases. For example, a staggered 
AB sublattice potential breaking the inversion symmetry 
leads to a quantum valley-Hall effect [2-4 ; certain non- 
magnetic adatoms, i.e., indium and thallium atoms [5], 
can enhance the intrinsic spin-orbit coupling (SOC) in 
graphene to host the quantum spin-Hall state 0; some 
3(i or 5(i transition metal adatoms produce a quan- 
tum anomalous-Hall state [7-9 , a consequence of the 
interplay between the proximity-induced magnetization 
and extrinsic Rashba SOC [10 . It is noteworthy that 
the intervalley scattering in these studies is completely 
avoided by choosing appropriate adatom coverages with 
valleys being separated in the momentum space, e.g., one 
adatom in each 4x4 supercell of graphene (6.25% adatom 
coverage) . 

However, at some specific adatom coverages, i.e., one 
adatom in each 3n x 3n (n = 1,2,3...) supercell of 
graphene, valleys K and K' are folded into the T point. 
The resulting intervalley scattering becomes significant 
and drives the system into a trivial insulator ^1 . In 
realistic graphene samples, the precise control of period- 
ically distributed adatoms is impractical. The randomly 
distributed adatoms will inevitably cause mixture of dif- 
ferent coverages and introduce the intervalley scattering. 
Therefore, a crucial issue arises: taking the site random- 
ization of adatoms into account, is it possible to observe 
the topological phases experimentally? 

In this Letter, we show that the randomization of 
adatom distribution greatly weakens the intervalley scat- 
tering, but does not affect the induced SOCs and magne- 



tization, therefore increasing the possibility of realizing 
topological phases in graphene. By making use of a finite- 
size scaling method, we show that the randomization can 
induce a topological phase transition from a trivial in- 
sulator to a topologically nontrivial phase. Using the 
Landauer-Biittiker formula, we confirm our finding by 
computing the two-terminal conductance in the presence 
of periodically or randomly distributed adatoms. 

Quantum spin-Hall effect — We consider a graphene 
sheet adsorbed with diluted nonmagnetic atoms (e.g. 
indium or thallium), which prefer the hollow site of 
graphene [5 . For simplicity, we assume that the adatoms 
only interact with the surrounding six-nearest carbon 
atoms. Such an interaction not only enhances the in- 
trinsic SOC, but also generates an on-site potential (also 
known as the crystal field stabilization energy) on each 
of the influenced carbon atoms. This potential is a key 
factor to induce the intervalley scattering. 

In the presence of randomly distributed adatoms, the 
tight-binding Hamiltonian of graphene can be written 
as [51 Hill: 

{ij),a {{ij))en,a(3 
+ ^ E 4aC^a. (1) 

where c]^ creates an electron on site i with spin a, and t is 
the hopping energy between nearest neighbors. The last 
two terms represent respectively the intrinsic SOC Aso 
and on-site energy [/, which are applied on the influenced 
atomic sites denoted by 1Z. Sz is the z-component of spin 
Pauli matrices. ((...)) sums over all next-nearest neigh- 
bors. Uij = 1( — 1) corresponds to the hopping clockwise 
(counterclockwise) between next-nearest neighbors. Ac- 
cording to their formation mechanisms, it is known that 
Aso should be one order smaller than U. In the following, 
we adopt the strength of intrinsic SOC in the thallium- 
atom adsorption case, i.e., Aso = O.OWt ~ 0.044 eV [5], 
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FIG. 1: (Color online) (a)-(b) Two different configurations of 
4 adatoms in a 6 x 6 supercell. (c)-(d) Bulk band structure 
in the presence of only on-site energy U — 0.36t [panel (c)] 
or intrinsic SOC Aso = 0.016t [panel (d)]. (e) Z2 topological 
number versus U at fixed Aso = O.OlGt. In (c)-(e), solid and 
dashed lines correspond respectively to the supercells in (a) 
and (b). (f) Schematic of the finite-size scaling at a fixed 
11.1% adatom coverage. Empty, small solid and large solid 
circles represent the pristine graphene, carbons influenced by 
adatoms, and adatoms, respectively. 



and the on-site energy is set to be = 0.36t ^1.0 eV. 

As mentioned in the Introduction, the topological 
phases in graphene are sensitive to the coverage of pe- 
riodically distributed adatoms. To make our investiga- 
tion complete and convincing, we shall discuss the effect 
of randomization on two kinds of adatom coverage with 
and without intervalley scattering. In the following, us- 
ing a finite-size scaling method, we simulate a fully ran- 
dom adsorption system by randomizing the distribution 
of adatoms in a gradually increasing supercell at the fixed 
adatom coverage. We begin with the discussion on the 
11.1% adatom coverage with intervalley scattering. 

Let us first examine the sensitivity of topologically triv- 
ial/nontrivial phases to the adatom configuration. Fig- 
ures [TJa) and [TJb) display two different 6x6 supercells 
of graphene. To keep the coverage, four adatoms are in- 
cluded in each supercell. In Fig.[TJa) the 3x3 periodicity 
still holds, while in Fig. [ijb) the 3x3 translational sym- 
metry is broken due to the shift of one adatom to its 
neighboring site. This can be regarded as the simplest 
step to randomize a 3 x 3 supercell to a 6 x 6 one. To 
understand the role of randomization on the on-site po- 
tential U and intrinsic SOC Aso 7 we plot the bulk band 
structure around F point by considering only U or Aso- 
In the presence of only [/, the intervalley scattering in 



the 3x3 supercell opens a bulk gap as plotted in solid 
line in Fig. [TJc) . One can see that a slight change of the 
adatom configuration shrinks the bulk band gap, imply- 
ing the decrease of the invtervalley scatterng. As a sharp 
contrast, the band structures in the presence of only in- 
trinsic SOC are almost identical for both supercells [see 
Fig. [TJd)]. This indicates that the effect of intrinsic SOC 
is insensitive to the slight change of the adatom configu- 
ration. 

Although both on-site potential and intrinsic SOC 
open bulk gaps, the resulting insulators are topologically 
different. An efficient way to identify these phases is 
to calculate the Z2 topological number, which charac- 
terizes the bulk band topology of a time-reversal invari- 
ant system [13 . By computing Z2 topological number 
using the methods discussed in Ref. [14], we find that 
when both U and Aso are nonzero, their competition re- 
sults in different phases. As shown in Fig. [ije), for a 
fixed Aso = 0.016t, small U leads to a Z2 = 1 quan- 
tum spin-Hall insulator, while large U gives rise to a 
Z2 = trivial insulator. A remarkable difference be- 
tween the two configurations can be observed in the range 
of G [0.29t, 0.36t]: a phase transition occurs from a 
trivial insulator to a Z2 = 1 topological insulator. 

Through employing a finite-size scaling method, we 
study the effect of strong randomization of adsorption 
sites. The numerical simulation procedures can be sum- 
marized as: (1) Expanding the supercell from 3 x 3 to 
3n X 3n (n = 2,3,4...); (2) Randomly selecting hol- 
low sites and then determining the influenced atomic 
sites IZ] (3) Based on the Hamiltonian in Eq. [l] calcu- 
lating the bulk band structure and measuring the bulk 
gap A; (4) Repeating the procedures (2) and (3) to ob- 
tain M samples. After that, the probability distribu- 
tion of bulk gap P( A) can be obtained using the formula 
P(A)JA = m/M, where m counts the magnitudes locat- 
ing within the range of [A - ^ A/2, A -h (5 A/2). 

Figure [2]exhibits the evolution of the probability distri- 
bution of bulk gap P(A) along with the increasing of the 
supercell size n. The left column corresponds to the case 
with only the on-site potential U = 0.36t. For n = 1, 
although there are Cg = 9 different adatom configura- 
tions, their band structures are exactly the same. There- 
fore, the resulting bulk gap is a constant A = 0.029t [see 
Fig. [ij^a)]. When n > 1, there are Cg^2 adatom config- 
urations in a 3n X 3n supercell, e.g., a 9 x 9 supercell 
has C^^ = 2.61 X 10^^ adatom configurations, most of 
which result in distinct band structures. In Figs. [2|b)- 
|2|e), we observe that the band gap A fluctuates in a 
wide region, and the gap region shrinks towards zero 
for larger n. To better reflect this characteristic, we in- 
troduce the median of each ensemble A, which is high- 
lighted in a red line. The exponential decay of A as 
a function of n can be easily visualized. In the inset of 
Fig.[2|e), we show that the probability distribution P(A) 
of n = 9 can be fitted by a 2D Maxwell distribution func- 
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FIG. 2: (Color online) Probability distribution of bulk gap P(A) for 3n x 3n supercells of graphene subjected to randomly 
distributed adatoms. 2 000 samples are collected for each panel, (a)-(e) Only the crystal field stabilization energy U = 0.36t is 
considered, (f)-(j) Only the intrinsic SOC Aso = 0.016t is considered, (k)-(o) Both U = 0.36t and Aso = 0.016t are included. 
The red line labels the median of each ensemble A. Inset: Numerical fitting of the probability distribution for n = 9. 



tion /(A) = A/cr2exp(-AV2cr2), where a = A/V2 In 2. 
According to /(A), the probability of opening a band 
gap in the range of [3A, oo) is about 0.2%. Together 
with the fact that a realistic graphene sample resembles 
a n ^ oo supcell, one can reasonably conclude that the 
intervalley scattering should be vanishing in the presence 
of randomly distributed adsorbates. 

Next we turn to the case with only the intrinsic SOC 
Aso = 0.016t. As shown in the middle column, we find 
that the bulk gap A is almost independent of the ran- 
domization for any supercell size n. This means that the 
intrinsic SOC is still insensitive to the strong randomiza- 
tion of adatom distribution, and the resulting quantum 
spin-Hall phase is stable against the random adsorption. 
Comparing the results from both limits of only on-site 
potential or intrinsic SOC, it is natural to expect that 
the realistic graphene sample should favor the quantum 
spin-Hall state when both U and Aso are present. 

Such a speculation is verified by the gap statistics 
drawn in the right column. In the 3x3 supercell, be- 
cause of the competition between U and Aso, the gap 
opening is smaller than those shown in Figs. [2|a) and 
M^f). When the supercell size increases, the median A 
first decreases at n = 2, but then increases for larger n 
[see Figs. |2|m)-[2|o)]. Through computing the Z2 topo- 
logical number, we obtain that Z2 = for n = 1, while 
Z2 = 1 for n = 3, 6, and 9. Therefore, it is clear that 
the increasing of randomization turns graphene from a 
trivial insulator to a quantum spin-Hall insulator. From 
the tendency of A, one can notice that it gradually ap- 



proaches the band gap labeled in Fig. |2jf). 

Thus far, we analyze the influence of randomization 
by investigating the bulk band structures of various su- 
percells. Hereinbelow, we design a two-terminal device 
illustrated in Fig.|3]to study the transport properties us- 
ing the Landauer-Biittiker formula [15 . The adatoms are 
only considered in the central scattering regime, and the 
two leads are modeled by pristine graphene ribbons. The 
presence of helical edge modes is one of the most striking 
properties of the quantum spin-Hall effect, which gives 
rise to a quantized longitudinal conductance. 

Figure |3|c) plots the average conductance G and its 
fluctuation SG as a function of the Fermi energy sf. The 
parameters are U = 0.36t and Xso = 0.016t. In the pres- 
ence of periodically distributed adatoms [see Fig. [sj^a)], 
G = and SG = in units of e'^ /h within the range of 
€f ^ [0.1 17t, 0.124t], signaling a trivial insulator. How- 
ever, when the adatoms become randomly distributed in 
Fig. [sj^b), a quantized plateau G = 2e^/h with vanishing 
fluctuation emerges in the regime of G [0.116t, 0.132t]. 
The distribution of local currents illustrated in the inset 
of (c) further indicates a quantum spin-Hall insulator. 
Such a phase transition resembles the disorder-induced 
topological Anderson insulator [l6l[T7]- The major dif- 
ference is that it is the adatom configuration but not the 
disorder strength to trigger the phase transition. 

Throughout the above analysis, we focus only on the 
11.1% adatom coverage with strong intervalley scatter- 
ing. What happens for other coverages without inter- 
valley scattering? Let us take the 6.25% adatom cover- 
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FIG. 3: (Color online) (a)-(b) Schematic of a two-terminal setup in the presence of periodically and randomly distributed 
adatoms, respectively. And both adatom coverages are 11.1%. (c)-(e) Comparison of conductances G between periodic and 
random adsorption as a function of Fermi level £f- (c)-(d) The on-site potential and intrinsic SOC are set ioheU — 0.36t and 
Aso = 0.016t. And the adatom coverages are 11.1% in panel (c) and 6.25% in panel (d). In (e), the on-site potential, Rashba 
SOC, and exchange field are set to be [/ = 0.38t, \r = 0.04t and M = 0.04t, respectively. And the adatom coverage is 11.1%. 
Circle and triangle symbols represent the periodic and random adatoms. The associated error bar denote the conductance 
fluctuation 5G. Insets of (c) and (e): schematic of the quantized channels from left to right. 



age (one adatom in a 4 x 4 supercell) for example. Us- 
ing the same finite-size scaling method, we show that 
the bulk gap is only dependent on the intrinsic SOC, 
but independent of the on-site energy or supercell size. 
An immediate evidence is the robust quantized plateau 
shown in Fig. [sjd) for either periodically or randomly 
distributed adatoms. To conclude, in a thallium- atom 
adsorbed graphene, the quantum spin-Hall state is a 
system-preferred ground state for any adatom coverage, 
and the bulk gap ~ 46 meV can be detected under cur- 
rent experimental techniques. 

Quantum anomalous Hall effect — When the 3d/5d 
transition metal atoms are adsorbed on graphene, the 
interaction enhances the Rashba SOC A^^ [18^, and in- 
duces the on-site potential U and magnetization M. In 
Refs. [8 and [10 , it is found by ab initio calculations 
that quantum anomalous Hall phase can be produced in 
the 6.25% adatom coverage, but the trivial insulator is 
usually formed in the 11.1% adatom coverage due to the 
strong intervalley scattering. In Fig. [sj^e), we calculate 
the average conductance versus the Fermi energy ep at 
a 11.1% adatom coverage. The parameters are chosen to 
heU = 0.38t, Xr = 0.04t, and M = OMt. We observe 
that a = Oe'^/h in the range of ep e [0.125t, 0.127t] 
for the periodically distributed adatoms. However, for 
the randomly distributed adatoms, a quantized plateau 
G = 2e'^/h emerges in the range of £f e [0.127t, 0.137t]. 
The inset plots the schematic of the chiral edge modes 



from left to right, which indicates a quantum anomalous 
Hall insulator. This result further confirms our finding 
that the intervalley scattering is fragile, but the spin- 
orbit couplings and exchange field are robust against the 
randomization of adatom distribution. 

Summary — In the presence of periodically distributed 
adatoms, the interaction induced spin-orbit couplings can 
produce various topological phases when the intervalley 
scattering is completely suppressed. At certain adatom 
coverage, the strong intervalley scattering plays a dom- 
inant role and turns graphene into a trivial insulator. 
Using both finite-size scaling method and transport cal- 
culation, we show that when the adatom distribution 
becomes random, the intervalley scattering is weakened, 
but other quantities (e.g. spin-orbit couplings, and ex- 
change field) are not affected. This finding points out 
that the topological states are graphene- favored ground 
states in the presence of randomly distributed adtoms. 
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